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Abstract 

The statistics of domain walls for ground states of the 2D Ising spin glass with +1 and -1 
bonds are studied for L x L square lattices with L < 48, and p = 0.5, where p is the fraction of 
negative bonds, using periodic and/or antiperiodic boundary conditions. When L is even, almost 
all domain walls have energy Ed w = or 4. When L is odd, most domain walls have E^ w = 2. The 
probability distribution of the entropy, Sd w , is found to depend strongly on E^ w . When E& m = 0, 
the probability distribution of \Sdw\ is approximately exponential. The variance of this distribution 
is proportional to L, in agreement with the results of Saul and Kardar. For E^w = k > the 
distribution of Sd w 1S n °t symmetric about zero. In these cases the variance still appears to be 
linear in L, but the average of Sd w grows faster than y/L. This suggests a one-parameter scaling 
form for the L-dependence of the distributions of Sdw f° r k > 0. 

PACS numbers: 75.10.Nr, 75.40.Mg, 75.60.Ch, 05.50.+q 
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I. INTRODUCTION 



There continues to be a controversy about the nature of the Ising spin glass. The 
Sherrington-Kirkpatrick model,— with its infinite-range interactions between the spins, is 
described by the Parisi replica-symmetry breaking mean-field theory.^* 3 - To understand mod- 
els with short-range interactions on finite-dimensional lattices, however, it is necessary to 
include the effects of interfaces, which do not exist in a well-defined way in an infinite- 
range model. The droplet model of Fisher and Huse,^*& which starts from the domain-wall 
renormalization group ideas of McMillan£&£ and Bray and Moore , WV 2 and studies the 
properties of interfaces, provides a very different viewpoint on the spin-glass phase. 

In two dimensions (2D), the spin-glass phase is not stable at finite temperature. Because 
of this, it is necessary to treat cases with continuous distributions of energies (CDE) and 
cases with quantized distributions of energies (QDE) separately . 11 ^ 3 

In three or more space dimensions, where a spin-glass phase is believed to occur at finite 
temperature T, the general framework of thermodynamics requires that the CDE and the 
QDE should be treated on the same footing. The way this comes about is that in these 
cases the typical domain wall energy increases as a positive power of the size of the lattice. 
Thus the quantization energy becomes a negligible fraction of the domain wall energy for 
large lattices. All bond distributions behave in a qualitatively similar way, except that the 
QDE have finite ground state entropies.—^ 

Amoruso, Hartmann, Hastings and Moore^ have recently proposed that in 2D there is a 
relation 

where ds is the fractal dimension of domain walls, and 9e is the exponent which characterizes 
the scaling of the domain wall energy with size. For the CDE case, the existing numerical 
estimates of ds and 9e satisfy Eqn. (1). However, it is unclear if Eqn. (1) should continue 
to be correct when the scaling exponent for spin correlations, rj, is not zero. For the QDE, 
the current estimate o 1 ^ 1 ^ find r\ rs 0.14. 

In three dimensions it is known from the droplet theory ,^' 12 that for the QDE, which 
have a positive entropy at T = 0, in the spin-glass phase 

d s = 29 s . (2) 
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6s is the exponent for the scaling of domain wall entropy with size. Thus, for the QDE, 
Eqn. (1) provides a relation between the scaling of the energy and the entropy of domain 
walls. It is not known how to calculate ds directly for the QDE case, so we need to use 
Eqn. (2) to check Eqn. (1) in that case. One might hope that this relation would also hold 
in 2D, even though the spin-glass order only occurs at T = 0. 

For the QDE, it is known that 6e = . 1Sl17 Then using Eqn. (1) gives ds = 5/4, or using 
Eqn. (2), 6 S = 5/8. The calculation of 6 S by Saul and Kardari&i 9 - found 6 S = 0.49 ± 0.02. 
Since ds cannot be less than 1, this result was interpreted as a strong indication that 
6 S = 1/2. 

In this work we will find that Eqn. (1) may not work for the QDE case in 2D. It appears, 
however, that Eqn. (2) is still correct in 2D, except when the domain wall energy, E^ w , 
is zero. The actual behavior of the QDE probability distributions under finite-size scaling 
turns out to be more subtle than what has been assumed until recentl y 20 i 21 As pointed 
out by Wang, Harrington and Preskill, 22 domain walls of zero energy which cross the entire 
sample play a special role when the energy is quantized. 

We will analyze data for the Edw and for the domain wall entropy, Sd w , for the ground 
states (GS) of 2D Ising spin glasses obtained using a slightly modified version of the computer 
program of Gallucio, Loebl and Vondrak,— which is based on the Pfaffian method. The 
Pfaffians are calculated using a fast exact integer arithmetic procedure, coded in C++. 
Thus, there is no roundoff error in the calculation until the double precision logarithm is 
taken to obtain Sdw- This extended precision is essential, in order to obtain meaningful 
results for entropy differences at large L. An earlier version of the domain wall entropy 
calculation,— using data provided by S. N. Coppersmith,— was limited to small L x L 
lattices with even L and came to somewhat different conclusions. 

We will demonstrate that for L x L square lattices the Edwards-Anderson 2 ^ (EA) model 
with a ±J bond distribution has a strong correlation between Ed w and Sdw for the GS 
domain walls. Because of this correlation, we will need to treat domain walls of different 
energies as distinct classes. We will find that the scaling parameter identified by Saul and 
Karda r 18 ' 19 is the one associated with domain walls having Ed w = 0. It is not, however, the 
one which controls the dominant behavior for large L. 
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The Hamiltonian of the EA model for Ising spins is 



H — — ^2 JijCTiCT 

(ij) 



3 ' 



(3) 



where each spin <jj is a dynamical variable which has two allowed states, +1 and -1. The {ij) 
indicates a sum over nearest neighbors on a simple square lattice of size L x L. We choose 
each bond Jjj to be an independent identically distributed quenched random variable, with 
the probability distribution 



so that we actually set J = 1, as usual. Thus p is the concentration of antiferromagnetic 
bonds, and (1 — p) is the concentration of ferromagnetic bonds. 

II. GROUND STATE DOMAIN WALLS 

We define the GS entropy to be the natural logarithm of the number of ground states. 
For each sample the GS energy and GS entropy were calculated for the four combinations of 
periodic (P) and antiperiodic (A) toroidal boundary conditions along each of the two axes 
of the square lattice. We will refer to these as PP, PA, AP and AA. In the spin-glass region 
of the phase diagram, the variation of the sample properties for changes of the boundary 
conditions is small compared to the variation between different samples of the same size,— 
except when p is close to the ferromagnetic phase boundary and the ferromagnetic correlation 
length becomes comparable to L. 

We define domain walls for the spin glass as it was done in the seminal work of McMillan.^ 
We look at differences between two samples with the same set of bonds, and the same 
boundary conditions in one direction, but different boundary conditions in the other direc- 
tion. Thus, for each set of bonds we obtain domain wall data from the four pairs (PP,PA), 
(PP,AP), (AA,PA) and (AA,AP). The reader should remember that the term "domain wall" , 
as used in this work, refers only to this procedure. Saul and Karda r 18 4 9 follow the same 
procedure used in this work, but use the term "defect" instead of "domain wall" . 

The domain-wall renormalization group of McMillan^ is based on the idea that we are 
studying an effective coupling constant which is changing with L. For the CDE case we can 
use the energy as the coupling constant. For the quantized energy case, what we need to 



P{J ij )=pS(J ij + l) + {l-p)5{J ij -l) 



(4) 
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do is a slight generalization of this idea. We should think of the coupling constant as the 
free energy at some infinitesimal temperature. When we do this, the entropy contributes to 
the coupling constant. As we will see, the distribution of E dw rapidly becomes essentially 
independent of L as L becomes large, except that there are separate distributions for even 
L and odd L. Under these conditions, it becomes possible to treat each value of E^ w as a 
separate class, representing a different coupling constant. 

The domain wall entropy, Sd w , is defined, by analogy to Ed w , to be the difference in the 
GS entropy when the boundary condition is changed along one direction from P to A (or 
vice versa), with the boundary condition in the other direction remaining fixed. [Sdw], where 
the brackets [ ] indicate an average over random samples of the J^-, is expected to increase 
as a positive power of L for any Ed w . Therefore, these coupling constants must eventually, 
at large enough L, be controlled by [Sdw] f° r an Y T > 0. Of course, the value of L which is 
needed for this to happen depends in T. The droplet model assumes that all these coupling 
constants, except for the Edw = case which has a special symmetry, are equal. 

As long as Edw > 0, the two boundary conditions which we are comparing are not on 
an equal footing. As Wang, Harrington and Preskill— express the situation, the Ed w > 
domain wall does not destroy the topological long-range order. However, in the Ed w = case 
the two boundary conditions are on an equal footing, and the topological order is destroyed. 
Therefore the Ed w = class of domain walls can be expected to behave in a special way, 
which differs from the prediction of the droplet model. 

It is natural to wonder if topological long-range order can be related to replica-symmetry 
breaking, and if the Edw = domain walls can be described by the replica-symmetry breaking 
theory. We will not attempt to do this here. 

It is important to realize that the meaning of a domain wall is very different when the 
GS entropy is positive, as in the model we study here, as compared to the standard case of a 
doubly degenerate ground state. In the standard case one can identify a line of bonds which 
forms a boundary between regions of spins belonging to the two different ground states. It 
is not possible, in general, to do that when there are many ground states. Despite this, we 
continue to use the term "domain wall" . 

When L is even, the energy difference, Ed w , for any pair must be a multiple of 4. When 
L is odd, Edw is 4n + 2, where n is an integer.— The sign of Ed w for a pair is essentially 
arbitrary for p = 1/2. Thus we can, without loss of generality, choose all of the domain- wall 
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TABLE I: Domain wall energy statistics for p = 0.5 with even L. The number of random bond 
configurations studied for each L was 500, and there are four McMillan pairs for each of these, rij 
is the number of domain walls of each type having Edw = k. fk is the fraction of domain walls 
having E dw = k. 



L n n 4 n 8 / / 4 

8 1467 530 3 0.7335 0.265 

12 1542 458 0.771 0.229 

16 1515 484 1 0.7575 0.242 

24 1578 422 0.789 0.211 

32 1530 470 0.765 0.235 

48 1546 450 4 0.773 0.225 



energies to be non-negative. 

III. NUMERICAL RESULTS 

Our calculated statistics for Ed w at p = 0.5, as a function of L, for even L and odd L are 
given in Table I and Table II, respectively. For each L, 500 distinct random configurations 
of bonds were studied. We obtain four McMillan pairs for each random sample, so we have 
2000 sets of Ed w and S dw at each L. For even L > 10 it turns out, crudely speaking, that 
about 77% of the time we find Edw = 0, and 23% of the time Ed w = 4. For odd L > 20, 
Edw — 2 about 98.5% of the time. No domain walls with energies greater than 8 were 
observed at any L for these values of p. This, however, does not have much fundamental 
significance. The probability distribution for Ed w is also a weak function of p— and a strong 
function of the aspect ratio of the lattice.— Our results are consistent with the results of 
Amoruso et alM 

It is interesting to note that Wang, Harrington and Preskill^ use an analytical argument 
to predict that f , the fraction of E dw = walls, is approximately 0.75, independent of p, 
in the spin-glass regime. However, the value of fo depends strongly on the aspect ratio of 
the lattice , 28 i 29 and it is not clear why this analytical argument should apply only when the 
aspect ratio is equal to one. It is also completely unclear to this author where the argument 
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TABLE II: Domain wall energy statistics for p = 0.5 with odd L. Column labels as in Table I. 



L n 2 n 6 f 2 

7 1944 56 0.972 

11 1960 40 0.980 

15 1957 43 0.9785 

21 1973 27 0.9865 

29 1967 33 0.9835 

41 1973 27 0.9865 




FIG. 1: (color online) (a) Average \Sd w \ vs. L for the E^ w = domain walls, log-log plot. The 
error bars indicate one standard deviation, (b) Histogram of \Sa w \ for E^ w = with L = 48. The 
vertical scale is logarithmic. 

uses the fact that E& w = domain walls can only occur when L is even. 

Estimating the statistical uncertainties in the data precisely is not trivial, due to the fact 
that the values of E^w obtained from the same set of bonds with the four different pairs of 
boundary conditions are not statistically independent M An upper bound on the statistical 
uncertainties is obtained by counting the number of samples, rather than the number of 
McMillan pairs of boundary conditions. 

The probability distribution of Sd w for the cases where Ed w = should be symmetric 
about 0, and our statistics are consistent with this. If we look at the L-dependence of 
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FIG. 2: (color online) (a) Average Sdw vs. L for the -E^™ > domain walls, log-log plot. The error 
bars indicate one standard deviation, (b) Variance of Sdw vs. L, log-log plot. 

shown in Fig. 1(a), we find a scaling exponent 

S (Q) = 0.500 ±0.020 (5) 

for Ect w = 0. The result of Saul and Kardar , 18 ' 19 obtained by looking at the distribution of 
Sdw f° r all values of Ed w combined, was 6$ = 0.49 ± 0.02. To obtain this exponent, Saul and 
Kardar fit their data at small values of Sdw- When L is even, which was the case for all of 
their data, this part of the data belongs almost entirely to the Ed w = component. 21 

The calculated means and skewness of these essentially symmetric distributions for Sdw 
is, naturally, consistent with zero, but their kurtosis is not. The reason for this is shown 
in Fig. 1(b), which is a histogram for \Sd w \ of Edw — when L = 48. We see that the 
distribution is approximately exponential, and therefore far from Gaussian. The computed 
kurtosis of this L = 48 distribution is 2.0, somewhat less than the value of 3 which would be 
found for an exact two-sided exponential distribution. The basic shape of these distributions 
is similar for the smaller values of L, with the width of each distribution given by the square 
root of its variance. 

When Edw is not zero, the relative signs of Ed w and Sd w are not arbitrary. Having chosen 
Edw t° be nonnegative, we then find that, when E dw is positive, it turns out that Sdw is 
usually positive. In Fig. 2(a) we show the behavior of [Sdw] for the cases where Edw = k, 
with k = 2, 4 and 6, as a function of L. We see that for k > 0, the average value of Sd w (L) 
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TABLE III: Scaling exponents for the first and second cumulants of the Sd w distributions. 9g is 
the scaling exponent for [Sdw], and (fts is the scaling exponent for the variance of (Sdw)- 



Edw 9s 4>s 

0.500 ± 0.020 0.992 ± 0.047 

2 0.565 ± 0.019 0.972 ± 0.051 

4 0.584 ±0.015 1.107 ±0.047 

6 0.617 ±0.062 0.85 ±0.28 



grows approximately as L . More precisely, least-squares fits to the form 

[S dw (L)\ = AL 9s (6) 

gives the results for 8s{k) shown in Table III. The result for k = 6 is rather uncertain, due to 
the small number of examples of this type. These results are consistent with the prediction 
of droplet theory^ that 8s should be independent of k (aside from the k = case, which is 
clearly exceptional). However, there also appears to be a tendency for 9s(k) to increase as 
k increases. Therefore, the possibility that 9s — > 5/8 as k — > oo, which would be consistent 
with Eqn. (1), cannot be excluded by these data. 

Because of the large GS degeneracy in the ± J Ising spin glass, one does not know how to 
compute ds directly for this model. However, if we use the droplet model prediction, that 
there is a single value for 9s, the result is not consistent with Eqn. (1). The author's opinion 
is that Eqn. (1) must be generalized when r\ > 0. 

As shown by Saul and Kardar , 18 ' 19 the variance of Sd w when Ed w = increases with L 
in approximately a linear fashion. Calculating the variance of these distributions, and using 
linear least squares fits on the log- log plot shown in Fig. 2(b), we find that assuming the 
increase of the variance with L is a power law gives the results shown in Table III. These 
numbers are reasonably consistent with the hypothesis that the scaling exponent for the 
variance of the Sdw distributions is equal to 1, independent of E dw . It is also interesting to 
observe that the magnitude of the variance, and not merely the slope of the fit, seems to be 
independent of E dw . Except in the special E dw = case, 29s is greater than 1. Therefore, 
the exponent ds should be controlled by 9s, as predicted by Eqn. (2). 
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FIG. 3: Histograms of S dw for (a) E dw = 2 with L = 41 and (b) E dw = 4 with L = 48. The 
vertical scales are logarithmic. 

IV. SCALING OF THE DISTRIBUTIONS 

In Fig. 3 we show histograms for the Sdw distributions for E dw = 2 at L = 41 and 
Edw — 4 at L = 48. In contrast to the E dw = case, the skewness and kurtosis of the 
Sd w distributions for E dw > are both small. It is possible that these distributions become 
Gaussian in the large L limit. However, the author is not aware of any reason why this must 
happen. 

The basic shapes of the histograms in Fig. 3(a) and Fig. 3(b) appear to be the same. 
Since 28 s > 1> it seems that the histogram for the E dw = 4 case can be mapped onto the 
histogram for the E dw = 2 case at a larger L. A way of expressing this is that for large L the 
Sdw histograms for E dw = k > should obey one-parameter scaling in the dimensionless 
variables 



[(Sdw) 2 ] - [Sdwl 



If 9s is independent of k, then, assuming <fts = 1 



g k {L) = {L/L k y e *-\ (8) 

where we define L k by the condition gk(Lk) = 1. 

What we have learned is that in this model there appear to be two distinct classes of 
domain walls, the E dw — domain walls and the E dv} > domain walls. As we have seen, 
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the Edw > domain walls behave in a way which appears to be essentially consistent with 
the predictions of the droplet model, but the E& w = domain walls do not. This difference 
in behavior is due to the symmetry of the E dw = case, which forces the average Sdw to be 
zero. 

For an E^ w > domain wall, a large contribution to Sd w comes from the shift in the 
average GS entropy with the shift in the GS energy— What remains to be understood is 
why [Sdw] should scale with L in the way predicted by the droplet model. The conventional 
derivation of the droplet models uses the assumption that the GS is unique, up to a reversal 
of the entire state, in an essential way. What follows immediately from this is that 77 = 
0. An extension of the droplet model to the more general case was given by Fisher and 
Huse. 6 However, the author hopes that by now he has convinced the reader that a better 
understanding of the r] > case is needed. 

V. SUMMARY 

We have studied the statistics of domain walls for ground states of the 2D Ising spin 
glass with +1 and -1 bonds for L x L square lattices with L < 48, and p = 0.5, where p is 
the fraction of negative bonds, using periodic and/or antiperiodic boundary conditions, for 
both even and odd L. Under these conditions, most domain walls have an energy Ed w < 8. 
The probability distribution of the entropy, S dw , is found to depend strongly on Ed w , but 
it appears possible to parameterize this dependence in a simple way. The results for Sd w 
do not appear to agree quantitatively with the prediction of Amoruso, Hartmann, Hastings 
and Moore,— Eqn. (1). Our results for [|Sd W |] when Ed w = agree with those of Saul and 
Kardar , 18 ! 19 but in addition we find that the distributions are close to being exponential in 
that case, even in the limit of large L. Due to the special role of the Edw — domain walls, 
we can understand the difference between the scaling exponent found by Saul and Kardar 
and the prediction of the droplet model. 
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